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Abstract 

The possibility of a strongly first-order electroweak phase transition is estab- 
lished in the minimal supersymmetric standard model with an extra U^l)', where a 
nontrivial CP violating phase is introduced in its Higgs sector. We find that there 
is a wide region in the parameter space of the model that allows the strongly first- 
order electroweak phase transition. The mass of stop quark need not be smaller 
than the top quark mass to ensure the first-order electroweak phase transition be 
strong. The effect of the CP violating phase upon the strength of the phase tran- 
sition is discovered. The strength of the phase transition is reduced when the size 
of the CP violation is increased. For a given CP violating phase, we find that the 
model has a larger mass for the lightest Higgs boson when it has a stronger phase 
transition. 
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I. INTRODUCTION 



The three necessary conditions for explaining the observed baryon asymmetry of the uni- 
verse, suggested several decades ago by Sakharov [1] , are the violation of baryon number 
conservation, the violation of both C and CP, and the deviation from thermal equilibrium. 
As is well known, in order to ensure sufficient deviation from thermal equilibrium, the 
electroweak phase transition (EWPT) that is responsible for the baryogenesis should be 
strongly first-order, since otherwise the baryon asymmetry generated during the EWPT 
would subsequently disappear. The criterion for strength of the EWPT is given by the 
baryon preservation condition Vc{Tc) > Tc [2], where Vc{Tc) and Tc are respectively the 
critical vacuum expectation value (VEV) and the critical temperature. 

In the literature, a large number of articles have been devoted to the possibility of 
establishing the strongly first-order EWPT in various theoretical models, including the 
Standard Model (SM) which is naturally the primary candidate. It has been found how- 
ever that the SM is inadequate to accommodate the desired EWPT, because it cannot 
produce the CP violation large enough to generate the baryon asymmetry by means of 
the complex phase in the Cabibbo-Kobayashi-Maskawa matrix alone [3,4], and it cannot 
make the EWPT strong enough unless the mass of the SM Higgs boson is below the 
present experimental lower bound. 

Attentions have been paid to alternative models beyond the SM within the context 
of the electroweak baryogenesis. Those models embedded in the supersymmetry (SUSY), 
having at least two Higgs doublets, are considered as the most plausible alternatives since 
they can bring about additional sources of CP violation in their Higgs sectors [5]. Among 
them, it has been shown that the minimal supersymmetric standard model (MSSM) 
with a light stop quark, being consistent with the present experimental bounds for both 
the Higgs bosons and the supersymmetric particles, allows the possibility of electroweak 
baryogenesis [6-14]. The tree-level Higgs potential of the non-minimal supersymmetric 
standard models, extended from the MSSM by introducing an additional Higgs singlet, 
may possess trilinear terms in Higgs fields. It has been studied that these trilinear terms 
may lead the models to the electroweak baryogenesis satisfying the baryon preservation 
condition in a wide parameter space, without requiring a light stop quark [15-21]. 

We have also studied several models for the possibility of the strongly first-order 
EWPT [22-24]. Recently, an MSSM with U{iy gauge symmetry has been our subject 
[24]. This model is found to satisfy the baryon preservation for a wide parameter space, 
without requiring scalar-pseudoscalar mixing or a light stop quark below top quark mass. 
This model may also accommodate new sources of CP violation in its SUSY soft terms, 
which allow enough baryon creation. It is mainly because its tree-level Higgs potential 
possesses the trilinear terms like the non-minimal supersymmetric models. 

In this paper, we continue to study the MSSM with f/(l)', this time including the 
explicit CP violation. In the explicit CP violation scenario, we examine carefully whether 
the strongly first-order EWPT may be realized in the model. We take into account the 
radiative corrections due to top and stop quarks, with the scalar-pseudoscalar mixings. 
We find that a strongly first-order EWPT is indeed viable for some parameter space of this 
model in explicit CP violation scenario. Then, we investigate the effects of the CP mixing 
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between scalar and pseudoscalar Higgs bosons, which is produced by the CP phase in the 
stop quark masses, on the strength of the first-order EWPT. We find that the CP phase 
arising from the stop quark masses reduces shghtly the strength of the first-order EWPT. 
On the other hand, the hghtest Higgs boson mass does not seriously affect the strength 
of the first-order EWPT, since it may be either enhanced or reduced by increasing the 
lightest Higgs boson mass, depending on other parameters. 



II. HIGGS SECTOR AT ZERO TEMPERATURE 

The Higgs sector of the MSSM with f/(l)' consists of two Higgs doublets Hi = {H^, H') 
and H2 = (H'^jH^), and a neutral Higgs singlet 5* [25-34]. Retaining only the Yukawa 
coupling of top quark, the superpotential of this model may be written as 

W ^ htQH2t% + XNH^eH^ , (1) 

where ht is the top Yukawa coupling coefficient, Q is the left-handed quark doublet su- 
perfield of the third generation, is the charge conjugate of the right-handed top quark 
superfield, A is the dimensionless coupling constant, and e is a 2 x 2 antisymmetric matrix 
defined as eu = 1. 

The tree-level Higgs potential of this model at zero temperature may be decomposed 
into F-terms, D-terms, and soft terms as 



Vo = VF + Vn + Vs, (2) 



with 
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Vf = \XWm'' + \H2\'')\S\^ + \H'[eH2 
Vd = 4iHlm + HlaH2f + 4i\Hi\'-\H2 



2\2 
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= mi\Hi\' + mi\H2\' + m'^\S\' -[XAxHleH2S + R.C.] , (3) 

where g2, gi, and g[ are respectively the gauge coupling coefficients of SU{2), U{1), and 
U{iy, a are Pauli matrices. Ax is the trilinear soft SUSY breaking parameter with mass 
dimension, rrii {i = 1-3) are soft SUSY breaking masses, and Qi [i = 1-3) are respectively 
the effective t/(l)' charges of Hi, H2, and S. 

The effective U (1)' charges of the Higgs fields satisfy the identity relation of Y^f^i Qi = 
in order to ensure the f/(l)' gauge invariance. The soft SUSY breaking masses are then 
eliminated by the three minimum conditions of the Higgs potential. The three minimum 
conditions define the vacuum and the vacuum expectation values of the three neutral 
Higgs fields at zero temperature, namely, vi = (Hi), V2 = {H2), and se*'^° = (S), where 
we assume that the Higgs singlet S may develop a complex VEV. 

Besides the phase 0^, the Higgs potential may explicitly contain complex coefficients. 
While other terms can be made real through redefinition of the Higgs fields, the trilinear 
terms, XAxHieH2S + B..C., may remain complex. Thus, we are left with a complex phase 
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00 = 0A + 0yiA where (px is the phase of A and is that of Ax- At the tree level, this 
complex phase 0o can always be eliminated by rotating suitably the relevant Higgs fields. 
The Higgs potential does not yield any CP mixing among the scalar and pseudoscalar 
Higgs fields [33,34]. 

At the one-loop level, a non-trivial complex phase emerges in the stop quark masses, 
which yields the mixing of the scalar and pseudoscalar Higgs bosons. In general, the 
contributions from the top and stop quark loops are most dominant at the one-loop level 
for a wide parameter space. The one-loop effective potential, at zero temperature, due to 
the top and stop quark loops is given as 



^i(O) 



^ 647r2 




(4) 



where A is the renormalization scale in the modified minimal subtraction scheme, and 
rif = —12 and n^. = 6 (z = 1, 2) are respectively the degrees of freedom for top and stop 
quarks from color, charge, and spin factors. 

We assume that the left-handed and right-handed stop quarks are not degenerate. 
Then, the stop quark masses are given as 



m 



tl,t2 



+ ^ J hfAfv2 + hfX^v(s'^ - 2hlXAtViV2S cos 



(5) 



where rriT is the soft SUSY breaking mass, rrit = htV2 is the mass of top quark. At is the 
trilinear SUSY breaking parameter with mass dimension, and 0t is a phase originating 
from the phases of A, At^ and s. 

In case of explicit CP violation, at the one-loop level, the non-trivial tadpole minimum 
condition with respect to the pseudoscalar component of the Higgs field is given as 



= A\ sin I 



3mfAt sin 0^ 



t2' 



(6) 



where 0o, which can be eliminated by rotating the neutral Higgs fields at the tree level, 
can not be so at the one-loop level, and the dimensionless function / is defined as 



m^, my) 



[mt 
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In the present model, there are four neutral Higgs fields. The physical Higgs bosons 
and their squared masses are given respectively by the eigenvectors and the eigenvalues, 
denoted as m\_ {i = 1-4), of the 4x4 mass matrix for the neutral Higgs fields. We sort 
these four neutral Higgs bosons in the increasing order of their masses such that m^^ is 
the smallest eigenvalue and hi is the lightest neutral Higgs boson. The explicit expression 
for the mass matrix may be found elsewhere [33,34]. 

The elements of the mass matrix, which are responsible for the mixing between the 
scalar and pseudoscalar Higgs bosons, are proportional to sin0t. Thus, the CP mixing in 
the present model at the one- loop level is triggered by 0^, which is practically the only 
source of CP violation in this model. If 0^ = 0, these four neutral Higgs bosons may be 
classified into three scalar and a pseudoscalar Higgs bosons. 
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III. THERMAL EFFECTS IN HIGGS POTENTIAL 

Up to now, we have studied the Higgs potential of the MSSM with f/(l)' at zero tem- 
perature. We now include the thermal effects to the Higgs potential in order to examine 
its behavior at the finite temperature. For the thermal effects, we take into account the 
one-loop contributions due to the top and stop quarks at finite temperature T, which is 
given as [36] 



1 + exp 



+ E 



j=l,2 



nrT^ 
27r2 



dx x^ log [l - exp (-y^x2 + A<|/r2j 



(7) 



where the first term represents the thermal effects of top quark, and second term those 
of the stop quarks, and the temperature dependence is explicitly shown. The masses of 
the field-dependent stop quarks in the above expression contain implicitly ipt that triggers 
the CP violation. The full effective potential at the one-loop level at finite temperature 
is therefore ^(r) = K) + ^1(0) + Vi{T). 

In order to study the possibility of a strongly first-order EWPT in this model, we 
examine if the full effective potential may develop two degenerate vacua at finite tem- 
perature, which is the critical temperature for the phase transition, and also examine if 
it may satisfy the baryon preservation condition. It is known that the trilinear terms 
proportional to A\ in the tree-level Higgs potential of the present model allows a strongly 
first-order EWPT for a wide parameter space, if the CP violation is turned off in its 
Higgs sector [34]. Then, we are mainly interested in the behavior of the present model 
in the presence of an explicit CP violation in its Higgs sector. In particular, we study 
the dependence of the strength of the strongly first-order EWPT upon the size of CP 
violation. 

In terms of the VEVs at finite temperature, the vacuum is determined by the minimum 
of the following effective potential. 



{V{v,,V2,s,T)) = {Vo) + (Vi(0)) + {V,{T)) , (8) 
where Vi, V2, and s are now the temperature-dependent vacuum expectation values with 



{Vo) 

{vm) 

{Vi{T)) 



,2„,2 



2„2 I fl'l + fi'i /„,2 



/ 2 2\2 I \2/' 2 2 I 2 2 I 2 2\ 



g ~ 

- 2XAxViV2S cos (j)o + -^(Qi-^i + Q2VI + Qss'^Y 



tt 327r2 
6T^ 



3m^ 
16^ 



TT^ Jo 



dx X log 



1 + exp (^-^Jx^ + ml{v2)/T' 



2 Q7^4 ^00 

+ E 



i=l 



dx x^ log 1 — exp + ''T^fX'^'i, V2, s)/T'^^ 



(9) 
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The minimum conditions for the one-loop effective potential at the zero temperature 
may express the soft SUSY masses mj (i = 1, 2, 3) in (Vq) as 

2 

mj = -^cos2/5-A2(s(O)2 + t;(O)2sin2/5)+AAAs(O)tan/5cos0o 
- g'lQiiQMOy cos^ f3 + Q2v{0f sin^ (3 + Q^siOf) 
+ -^{Xh{Of - \As{0) cot/5cos0,)/(m|,m|) , 



2 

2 

2~ 



ml = ^cos2/3 - X'^{s{Of + v{Oycos'^ 13) + XAxs{0) cot p cos (f)o 



- 9iQ2{QM0y cos' P + Q2v{0y sin' (3 + Q^siOY) 

+ _ AA,s(0) tan/3cos0i)/K ) - log ( ^^f^ 



= -A2i;(O)' + ^^t;(O)MAsin2/5cos0o 

- g'lQsiQMOy cos' (3 + g2t^(0)' sin' /3 + 03^(0)') 

where Wi(0), ^2(0), and s(0) are the zero-temperature vacuum expectation values of the 
three neutral Higgs fields. Notice the presence of the phase 0o; which is expressible as the 
other parameters via the tadpole minimum condition at the one-loop level shown in the 
previous section. 

The procedure of exploring the parameter space of the present model starts with 
reducing the number of variables. By calculating dV/ds\y-^^y^^s = 0, which is one of the 
minimum conditions that the first derivatives of the effective potential should satisfy, we 
obtain the following nonlinear equation, 

i2/„.2 



= 2m^^s - 2XAxViV2 cos 0o + 2A'K + ^2)^ + 2g/Q3s{Qivf + Q2V2 + Q3S 



^^^^{Xvis - AtV2 cos 0t)/(mf^,m|) 



3T^ 2hfXvi 
2ti'^ im? — 



— mf ^ 



exp(-w'x2 + m?yT2) 

X dx X 



x^ + ml/T^{l - exp(- + m|/T2)} 



3T' 2h\Xv^ ,^ 

+ TTT^^ ttKASVx - AtV2 cos (pt) 

27r2 (mf — mj ) 



/•oo exp(— ,/x2 + m? /T2) 
X dx x' . ^ , (111 



x2 + m? /T2{1 - exp(-./x2 + m? /T^)} ' 



which can be used to express s in terms of other variables, Vi and V2. Thus, by substituting 
s into {V{vi,V2,s,T)), we are left with (V(t>i, t>2, T)). 
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For given set of relevant parameter values, we search the minima of {V{vi,V2,T)) in 
the (ui, W2)-plane, varying the temperature. If two independent but degenerate minima 
are found for a given temperature, it is the critical temperature, T^, and the shape of the 
potential at the temperature allows a first-order EWPT. Let one minimum point in the 
(fi, f2)-plane be A, and the other, B. The distance between A and B is given as 

Vc = \J {VIB - ViaY + {V2B - V2aY + {sb - SaY , (12) 

where sa and s b are obtained by inverting the above nonlinear equation, by means of the 
bisection method. The critical vacuum expectation value, fc, is defined as the distance 
between the two degenerate minima. Then, by comparing Vc with the critical temperature, 
we can check if the baryon preserving condition is satisfied. In this way, the possibility 
of a strongly first-order EWPT is established in the MSSM with U{1)\ with explicit CP 
violation. 

IV. NUMERICAL ANALYSIS 

Now, for numerical analysis, we set the renormalization scale A as 300 GeV and the mass 
of top quark as 175 GeV. The effective U{1)' charges, Qj, may be redefined as Qi = giQi, 
for the sake of convenience, because Qi always go together with the U{iy gauge coupling 
constant [34]. Like Qi, Qi also satisfy the 11(1)' gauge invariance condition, J2Qi = 0- We 
have established the experimentally allowed region in the (Qi, (52)-plane in Ref. [34], using 
such constraints as the extra gauge boson mass and the mixing angle between the two 
neutral gauge boson azz'- We adopt the results of Ref. [34], namely, Qi = —I, Q2 = —0.1, 
and Q3 = 1.1. We further set for our numerical analysis tan/5 = 3, s(0) = rriT = 500 
GeV, and At = 100 GeV. Then, we are left with four free parameters: (pt, A, Ax, and T 
in the full effective potential. 

In Fig. 1, the equipotential contours of the effective potential at the one- loop level, 
{V{vi,V2,T)), are plotted in the (ui, Va) -plane for (pt = tt/IOOO, A = 0.8, Ax = 2271 GeV, 
and T = 100 GeV. We take a very small value for the phase (f)t in order to study the 
case of an almost exact CP conservation. Since there are two degenerate minima in Fig. 
1, T = 100 GeV is the critical temperature, and a first-order EWPT is possible. Let us 
examine the strength of the phase transition. The coordinates of the two minima in the 
(vi, V2)-plane determine the VEVs at the two degenerate vacua as: (w 1, V2, s) = (1, 3, 489) 
and (199, 507, 617) in GeV. Note that s is obtained indirectly through the nonlinear 
equation. The critical VEV is then calculated to be about Vc = 557 GeV for Fig. 1. 
Thus, the baryon preserving condition is satisfied enough since Vc/Tc = 5.57. It is therefore 
possible that the present model accommodates a strongly first-order EWPT. The masses 
of the four neutral Higgs bosons are obtained using the parameter values in Fig. 1 as 73, 
792, 1748, and 1750 GeV. 

In Fig. 2, the equipotential contours of {V{vi,V2,T)) are plotted in the (f 1, t;2)-plane 
for 0t = it/2, a = 0.8, Ax = 2271 GeV, and T = 147 GeV. We take the maximal value 
of the phase (pt for the CP violation, that is, sin0j = 1, in order to study the effect of 
the CP phase on the strength of the phase transition. The critical temperature in Fig. 2 
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is 147 GeV, since there are two degenerate minima. The coordinates of the minima are 
the same as Fig. 1, hence the critical VEV is the same as Fig. 1. However, the critical 
temperature is different from Fig. 1, and therefore the strength of the phase transition is 
different from Fig. 1: it is Vc/Tc = 3.79. We find that the strength of the phase transition 
is reduced. Also, we find that the mass of the lightest neutral Higgs boson is slightly 
increased as (pt is changed from tt/IOOO to 7r/2. 

We continue our numerical analysis for other parameter values. In Fig. 3, the equipo- 
tential contours of {V{vi, V2, T)) are plotted in the (wi, V2)-pla.ne for (pt = vr/1000, A = 0.7, 
Ax = 2115 GeV, and T = 100 GeV. The strength of the first-order EWPT is estimated to 
be about 6.24. The mass of the lightest neutral Higgs boson is about 79 GeV. In Fig. 4, 
for (f)t = tt/2, a = 0.7, Ax = 2115 GeV, the critical temperature is about Tc = 143 GeV. 
We find that the strength of the first-order EWPT is decreased from 6.24 to 4.36, and 
the mass of the lightest neutral Higgs boson is slightly increased from 79 GeV to about 
80 GeV, as (pt is changed from vr/lOOO to 7r/2, comparing the numerical results of Fig. 4 
with those of Fig. 3. 

The effect of (pt that accounts for the CP mixing between the scalar and pseudoscalar 
Higgs bosons upon the strength of the first-order EWPT seems rather straightforward, 
if we compare Fig. 1 with Fig. 2, or Fig. 3 with Fig. 4. Increasing the size of the 
CP mixing between the scalar and pseudoscalar Higgs bosons would result in decreasing 
the strength of the phase transition. On the other hand, the relationship between the 
strength of the first-order EWPT and the mass of the lightest Higgs boson may be seen 
clearly if we compare Fig. 1 with Fig. 3 or Fig. 2 with Fig. 4, since in this way we can 
exclude the effect of (pf Comparing Fig. 1 with Fig. 3, where (pt is fixed at vr/lOOO, it 
is evident that the strength of the phase transition increases as the mass of the lightest 
Higgs boson increases. This behavior is also observed if we compare Fig. 2 with Fig. 4, 
where (pt is fixed at 7r/2. Combining these two observations, one might expect that in the 
MSSM with U{iy in explicit CP violation scenario the strength of the phase transition 
increases if the mass of the lightest neutral Higgs boson increases or if the size of the CP 
mixing between the scalar and pseudoscalar Higgs bosons decreases. We summarize in 
Table 1 the relevant numbers for Figs. 1-4. 

V. CONCLUSIONS 

In the MSSM with U{1)' in explicit CP violation scenario, at the effective one-loop level, 
we study the possibility of a strongly first-order electroweak phase transition. The radia- 
tive corrections are evaluated by employing the effective potential method, taking into 
account the top and stop quark loops. The explicit CP violation is triggered by a phase 
arising in the stop quark masses. We have searched for the parameter space of the model 
to find that some region in the parameter space allows a strongly first-order electroweak 
phase transition. Four sets of parameter values are selected and plotted in figures in or- 
der to show explicitly that a strongly first-order electroweak phase transition is possible 
indeed. 

It is natural that the strength of the phase transition would change if we choose 
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TABLE 1: The values of relevant parameters in Figs. 1-4 which allow a strongly first- 
order EWPT in the MSSM with t/(l)' in explicit CP violation scenario, and the numerical 
results. The remaining parameters are set as tan/3 = 3, s(0) = mj- = 500 GeV, and 
At = 100 GeV. The arrow in a cell indicates that the values in the cell is equal to those 
in the adjoining left cell. 
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iviA,V2A,SA) (GeV) 


(1,3,483) 




(2,4,484) 




{vib,V2b,sb) (GeV) 


(199,507,617) 




(230,564,639) 




Vc/Tc 


5.57 


3.79 


6.24 


4.36 



different sets of parameter values. In particular, it may depend on the size of the phase 
for the CP mixing between the scalar and pseudoscalar Higgs bosons. We find that the 
strength of the strongly first-order electroweak phase transition decreases as the phase for 
the CP mixing increases. 

The mass of the lightest neutral Higgs boson is also changed as the parameter values 
are changed. We find that, for a fixed value of the CP phase, the lightest neutral Higgs 
boson becomes heavier when the strength of the phase transition is enhanced. In other 
words, if the phase transition is stronger for a set of parameter values than for another 
set of parameter values, then the mass of the lightest Higgs boson is also larger for the 
former set of parameter values than for the latter. 
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FIGURE CAPTION 

FIG. 1. : Equipotential contours of {V{vi,V2,T)) in the (wi, U2)-plane for (pt = vr/lOOO, 
tan/3 = 3, A = 0.8, Ax = 2271 GeV, s(0) = tut = 500 GeV, and At = 100 GeV, at 
the critical temperature T = 100 GeV. The two minima in the {vi,V2)-pla.ne determine 
the VEVs at the two degenerate vacua as: {vi,V2,s) = (1,3,489) and (199,507,617) in 
GeV. The strength of the first-order EWPT is about Vc/T^ = 5.57. The masses of the 
four Higgs bosons are nih^ = 73 GeV, rrih^ = 792 GeV, rrih^ = 1748 GeV, and m/^^ = 1750 
GeV. 

FIG. 2. : Equipotential contours of {V{vi,V2,T)) in the (f i, f2)-plane for (pt = vr/2, 
tan/3 = 3, A = 0.8, Ax = 2271 GeV, s(0) = = 500 GeV, and At = 100 GeV, at the 
critical temperature Tc = 147 GeV. The VEVs at the two degenerate vacua are the same 
as Fig. 1. The strength of the first-order EWPT is about Vc/Tc = 3.79. The masses of the 
four Higgs bosons are ruh-^ = 75 GeV, = 792 GeV, ruhs — 1747 GeV, and m/i^ = 1749 
GeV. 

FIG. 3. : Equipotential contours of {V{vi,V2,T)) in the (wi, f2)-plane for (pt = tt/IOOO, 
tan/3 = 3, A = 0.7, Ax = 2115 GeV, s(0) = tut = 500 GeV, and At = 100 GeV, at 
the critical temperature = 100 GeV. The two minima in the (t'l, f2)-plane determine 
the VEVs at the two degenerate vacua as: {vi,V2,s) = (2,4,484) and (230,564,639) in 
GeV. The strength of the first-order EWPT is about Vc/T^ = 6.24. The masses of the 
four Higgs bosons are m^^ = 79 GeV, = 789 GeV, m^g = 1579 GeV, and m/^^ = 1580 
GeV. 

FIG. 4. : Equipotential contours of {V{vi,V2,T)) in the (t^i, t'2)-plane for (pt = vr/2, 
tan/3 = 3, A = 0.7, Ax = 2115 GeV, s(0) = tut = 500 GeV, and At = 100 GeV, at 
the critical temperature Tc = 143 GeV. The VEVs at the two degenerate vacua are the 
same as Fig. 3. The critical temperature is determined by Tc = 143 GeV. The strength 
of the first-order EWPT is about Vc/Tc = 4.36. The masses of the four Higgs bosons are 
ruh, = 80 GeV, rrih^ = 789 GeV, mh., = 1578 GeV, and m/,, = 1579 GeV. 
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FIG. 1: Equipotential contours of {V{vi,V2,T)) in the (f i, W2)-plane for 0t = tt/IOOO, 
tan/3 = 3, A = 0.8, Ax = 2271 GeV, s(0) = ttit = 500 GeV, and At = 100 GeV, at 
the critical temperature T = 100 GeV. The two minima in the (wi, W2)-plane determine 
the VEVs at the two degenerate vacua as: {vi,V2,s) = (1,3,489) and (199,507,617) in 
GeV. The strength of the first-order EWPT is about Vc/Tc = 5.57. The masses of the 
four Higgs bosons are m^^ = 73 GeV, nih^ = 792 GeV, m^g = 1748 GeV, and nih^^ = 1750 
GeV. 
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FIG. 2: Equipotential contours of {V{vi,V2, T)) in the (f i, f2)-plane for (pt = 7r/2, tan (3 = 
3, A = 0.8, Ax = 2271 GeV, s(0) = ttit = 500 GeV, and At = 100 GeV, at the critical 
temperature Tc = 147 GeV. The VEVs at the two degenerate vacua are the same as Fig. 
1. The strength of the first-order EWPT is about Vc/Tc = 3.79. The masses of the four 
Higgs bosons are nih^ = 75 GeV, nih^ = 792 GeV, rrih^ = 1747 GeV, and m^^ = 1749 
GeV. 
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FIG. 3: Equipotential contours of {V{vi,V2,T)) in the (f i, W2)-plane for 0t = tt/IOOO, 
tan/3 = 3, A = 0.7, Ax = 2115 GeV, s(0) = ttit = 500 GeV, and At = 100 GeV, at 
the critical temperature Tc = 100 GeV. The two minima in the (f i, f2)-plane determine 
the VEVs at the two degenerate vacua as: {vi,V2,s) = (2,4,484) and (230,564,639) in 
GeV. The strength of the first-order EWPT is about Vc/Tc = 6.24. The masses of the 
four Higgs bosons are m^^ = 79 GeV, nih^ = 789 GeV, nih^ = 1579 GeV, and m^^ = 1580 
GeV. 
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FIG. 4: Equipotential contours of {V{vi,V2, T)) in the (f i, U2)-plane for 0^ = n/2, tan P = 
3, A = 0.7, Ax = 2115 GeV, s(0) = mr = 500 GeV, and At = 100 GeV, at the critical 
temperature Tc = 143 GeV. The VEVs at the two degenerate vacua are the same as Fig. 
3. The critical temperature is determined by Tc = 143 GeV. The strength of the first- 
order EWPT is about Vc/Tc = 4.36. The masses of the four Higgs bosons are rrih^ = 80 
GeV, nih^ = 789 GeV, = 1578 GeV, and m^^ = 1579 GeV. 
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